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I.  Introduction 

The  elm  of  the  present  program  was  to  program  the 
Stochastic  Process  Method  (SPIO^  for  neutron  problems  and  to 
Investigate  the  method  both  theoretically  and  numerically.  We 
have  been  able  to  analyze  the  problem  sufficiently  to  determine 
proper  source  and  Interface  conditions,  as  well  as  appropriate 
methods  of  treating  scattering  by  hydrogen  and  Inelastic 
scattering  by  heavy  elements  within  the  spirit  of  the  method. 

In  this  report  are  presented  the  SPM  equations  with 
boundary  conditions,  a  discussion  of  some  of  the  results,  and 
operating  instructions  for  the  IBM-7090  SPM  code  and  for  the 
auxiliary  Input  code  GENUC. 


Manuscript  released  by  the  authors  April  1962  for  publication  as  an  ASD 
Technical  Documentary  Report. 
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II.  Analysis  of  the  SPM  Equation 


A.  Primitive  SPM  Equation  for  Variable  Cross-Sections 

An  SPM  equation  can  be  derived  for  cross  sections  which  vary 
in  space  and  energy.  The  previous  derivation  was  valid  only  for 
constant  cross  sections.  We  now  derive  the  new  result. 

The  integral  Boltzmann  equation  can  be  written 


a(x,u)0(x,u)  =  1^q(x,u)  +  2tf 


u  rl 

fldxA  du' \  dn  I  dp,* 

-'-1  J-1 


g(x.u)  „/X-X'v 
ini  n  } 


X  exp 


/(■''.  "I 

I  -  \  a(x”,u)dx”  «•  \  dv\  d0  6  ju-u’-L(v)J  x 

X  6(p.-M.'v  -  ^1-p.*^  ^l-v^cos0)  Og(x»,uW)0'’(x‘,u»,p.*) .  (i) 


Here 


X  »  distance  from  source  plane 
u  =  lethargy 

p.  =  cosine  of  angle  with  normal  to  source  plane 
V  =  cosine  of  scattering  angle 
0^=  azimuthal  angle  of  scattering 

o(x, u)  =  total  macroscopic  cross  section  at  x  and  u 

0  (x,u,v)  -  macroscopic  differential  elastic  scattering 

s 

cross  section  at  x  and  u  for  scattering  through  an  angle 
whose  cosine  is  v,  per  unit  solid  angle 
0(x,u,M.)  -  flux  at  (x,u,p.)  per  unit  lethargy  and  solid  angle 
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L(v)  ■  increase  in  lethargy  accompanying  scattering  through 
an  angle  whose  cosine  is  v 

^q(x,u)  =  density  of  first  collisions  at  x  and  u  per  unit 
volume,  time,  and  lethargy 

Primed  variables  refer  to  the  situation  before  collision. 
Also,  0(x,u)  is  the  total  flux  at  x  and  u  per  unit  lethargy: 

0(x,u)  =  21^  0(x,u,u)  dn  (2) 

and  H(^)  is  a  step  function: 

H(0  =  i  >  0 

=  0,  i  <  0.  (3) 

Let  the  total  number  of  mean  free  paths  at  lethargy  u  between 
0  and  X  be 


s(x,  u)  =  \  a(x",u)du". 

o 


(4) 


s  is  thus  a  function  of  x  with  parameter  u  and  conversely,  x  is  a 
function  of  s.  We  will  also  use  the  notation 


s'  «  s(x' ,  u)  . 


(5) 


Let 


a  (x'  ,u' ,v) 

-  0(x',u',M.')  =  f(x'(s'),  u',u,v,M.')  =»  g(s'). 

a  (x'  ,  u) 


(6) 
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By  Taylor’s  expansion^ 


g(s')  *  2  —  (s'-s) 

n=0  n! 


g(s). 


(7) 


where  by  (6) , 


a  (x,u',v) 

g(s)  »  f(x(s),  u»,u,v,>i»)  =  g  - 0(x,u’,ti’).  (8) 


We  have  also 


H(^)  =  H(^). 


Then  Che  x'  -  Integral  in  (1)  can  be  written 


(9) 


-  a(x'*,u)dx'‘|  ffg(x»,u,v)  0(x',u’,n’)dx*  = 

'x»  ^ 


g(x,u)  { dx-  H(^)  e  °o''(x»-;ar 


joo 

,f§r 


n  <ti 


—  Co 


(s’“S)  g(s) 

OS 


Z  (-l)"n”a(x,u)  7-rg(s)* 
n=0  ds 


(10) 


The  second  equality  in  (10)  holds  by  virtue  of  (4) ,  (6) ,  (7) , 
and  (9) .  But  by  (4) , 

so 

y  0(x,usn>)  - 

-[ir  (12) 
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Inserting  (12)  into  (10)  and  then  (10)  into  (1) ,  we  have 

-,i  ^2'Tr  ri 

,  ^v.iiV  4.  9'rr“  \  rfii»  \ 

'^O 


(x,u)0(x,u)  =  ^^^(x,u)  +21^2  (-1)"  1  du*  )  dv\  d0^\  dn»\  du»  x 

naO 


'-1 


^-1  -  0 


-I  n 


^  a(x,u) 


n 


B(u- 


“1  ^  0 

M-'v  -  ^l-n*^  ^ l-v^cos0)  6  |u-u»-L(v^  X 


X  Cg(x,u*,v)0(x,u»,M.')  •  (13) 

We  now  expand  00  about  u: 

^  ,  m  -vm  o  (x,u,v) 

ag(x,u»,v)0(x,u*,M.»)  =  77  (“'"“)  'g  (x/u)^  Og(x,u)0(x,u)M(M.*  |x,i. 


maO  m! 


au*"  °s' 


where  o. (x,u)  is  the  elastic  scattering  cross  section  at  x  and  u; 
s 


Og(x,u)  a  2^  Og(x,u,v)dv 


(14) 


(15) 


and 


Mfc-  |x.u)  .  .  (16) 

Inserting  (14)  into  (13),  integrating  over  u’  and  p.,  and  then 
dropping  the  prime  on  p,',  we  have  the  SPM  equation: 


o(x,u)  a  V'o(x.u)  + 


“  “  (-1)"^"  fa  1  ~l  ” 

nfo  m! 


Nm 


(17) 
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where 


=  2'n'\  M(nlx,u)dM.  \  l"*(v) 

-1  -^-1 


Cfg  (x,u,v) 
ag(x,u) 


ZTT' 


dv 


^  0 


J^v  l-v^cos0  !  d0. 


(18) 


Equations  (17)  and  (18)  are  precisely  the  same  relations  as  for  constant 
cross  sections,  except  that  now  we  must  make  sure  that  the  a(x,u)  are 
placed  as  given  in  (17),  since  a  no  longer  commutes  with  the  derivatives. 
If  the  medium  consists  of  a  number  of  nuclides,  each  has  its  own 
The  appropriate  generalization  retains  eq.  (17)  with  defined 
by  the  relation 

(19) 


Ca  =2C^a^, 
mn  s  ^  mn  s' 


where 


1  ^  V 

a-  =  macroscopic  scattering  cross  section  for  i~  nuclide 

O 

‘^s  =  I 


and  is  defined  by  eq.  (19)  with  the  cross  sections  oj"  being  used. 

=  1. 


mn 

Note  that  C 


oo 

B.  The  Truncated  SPM  Equation  and  the  Expression  for  the  Current 

To  second  order  in  both  x  and  u,  the  SPM  equation  (17)  reduces 
in  the  absence  of  a  source  V'q  bo 

2 

.2 


35E  a  ^  ^11*^8^  ■  -55^  a  ^  a  ^  ^12  ®s^  "  7  35^  a  ^  ^21®s 


Id  Id  1  d^ 


7  ^  a  ''22  °s 


0^0  =  0 


(20) 


I 


It  was  shown  in  the  original  summary  report'*’  that  the  current  J 
obeys  the  equation 


hj  _  % 


Co 

z 

m=0  n=l 


llll 


m+n 


m: 


a  1 

^'Sx  a(x,  u) 


—  C  a  0, 
Su"  ® 


where  we  have  used  the  general  form  (17)  of  the  SPM  equation. 
In  the  order  to  which  we  are  truncating,  we  can  then  write 

s2 


1  r 

0  L 


^01  ^02  +  Is  ^11  •  llE  In  ^12  "  ■?  f;T  ^21 

I  h  I  a" 


.id  id  r  'oi 

■*'7  'S^a  ^  ^22  °s^* 


(22) 


Now  let 


(23) 


is  defined  here  somewhat  differently  from  the  way  in  which  it  was 
defined  previously.^ 

Then 


-  (1-v^,)  J  -  7  (v^lO  J)  +  75  (v2<yJ)  = 


I  a  1 

a  Sx  o 


=02  “si*  -  Is  <=12  “s®  +  1  It  <=22'>s« 


(24) 


Suppose  that  successively  higher  u-derivatives  correspond  to 
successively  higher  powers  of  some  small  parameter  "h.  Thus  if  some 
quantity  f  is  of  order  unity  [o(l)J^ 

B  ■ 

Expand  J  in  orders  of  ^  such  that  =  0(a  ) .  Then  to  second  order, 

J  »  Jq  +  +  ^2- 
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Again  to  second  order,  (24)  becomesi 


o  ^  a  I 


a 


1 


^02° a^  "  I  ^  (^22*^8^^ 


']• 


Then 


-  (l-V*^o  VF  lx  ?  ^02  ‘^s^ 

-  (l-v^)Ji  -  i  (v^a  J^)  =  -  i  i  C^2 


(26) 

(27) 

(28) 


(1-Vq)J2-  J  Ij  (v^ojp  +  ^  ^(V20J^)  =  J  i  ^  ^22(1^0. 


From  eqs.  (27) -(29),  we  get 

T  -  1  IS  1  r  ^  a 

•^o  ■  T=7:  0  ^  a 

Ji  * 


(29) 


(30) 


1 


1  [|>.  i  I-  Cl 0  +  I-  T-^  I-  “  0  (31) 

a  ^ox  a  ou  12  a^  ^  l-v^^  ox  a  02  s^J  ^  ' 


•^2  “  ■  T^  F  ^  li?  ^22®s^  I 

1 

^Fx  o 


C  0  0  + 
57  TTF;  Bx  o  Co2'^s'®  + 


BF  t::^ 


If  Ci2'^8^  +  If  T::^  I5E  F  ^02*^"^ 


■}] 


(32) 


o  \-  o 

Only  for  constant  cross  sections  can  expressions  (30) -(32)  be  com¬ 
bined  to  give  a  current  of  the  form 


In  the  general  case,  we  must  use  eqs.  (25),  (30),  (31),  and  (32)  to 
give  the  expression  for  J.  J  is,  of  course,  taken  to  be  continuous. 
We  obtain  the  truncated  SPM  equation  in  the  form 

2 

■  ■  lu  “l®^  “2^  ■  lx  “  ° 

where 

Oq  =  <T  -  Os  C35a) 

“l  “  ^10‘^s  (35b) 

“2  "  7  ^20*^8 '  (35c) 

Because  the  J2  term  was  found  not  to  improve  the  results,  the 
approximation  was  actually  cut  off  after 
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C.  Interface  Conditions 

The  conditions  at  an  Interface  must  be  continuity  of  0  and 
J,  for  all  lethargies.  It  is  necessary  to  see  how  this  works  out 
in  a  difference  approximation. 

We  assume  an  interface  at  x  =  0.  Then 

2 

"  “o^  "  lu  “l^  “2^  ‘  lx  “  ° 

We  will  denote  the  value  of  J  at  the  interface  by  J  .  If  6  .  is 

the  x-interval  to  the  first  mesh  point  to  the  right  of  the  interface 
and  6_  that  to  the  mesh  point  to  the  left,  we  have 

,  (37) 

where 

J_  =  value  of  J  at  the  first  mesh  point  to  the  left 
=  average  value  of  in  left-hand  Interval. 

The  mean  value  of  J  in  the  left-hand  Interval  is 

T:  (38) 

Then 

-  j:  =  —5^  -|a.  -  v  (39) 

Similarly, 
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If  is  the  value  of  to  the  right  of  the  Interface  and 
that  to  the  left,  we  have 

.2 


a 

a 


an 
.2 


■®+“o4^"  ^  ^  6^02/  +  2(J^  -  J^)  -  0 


(41a) 

(41b) 


Defining 


Oa 


6+0^++  6_aj,. 
6+  +  6. 


(42) 


and  adding  the  two  equations  (41  a,b)  we  have 

"“o^‘  In  “1^  ^  “2^  B"  l~bl  ■  ^4^ 


(43) 


Equation  (43)  is  precisely  what  we  have  within  a  single  medium. 
Thus  we  can  treat  boundary  points  the  same  as  any  other  points  if  we 
define  the  at  the  boundary  by  eq.  (42),  The  cross  section- 
dependent  constants  in  are  those  appropriate  to  the  right-hand 
medium  and  those  in  7.  are  those  appropriate  to  the  left-hand 
medium.  We  can  write 


-J'. 


\h  “3^  "  lx  ^4  lu  “4^  “  In  ^5  lx  “3^1  ' 


(44) 


where 


^3  “ 

1 

^4*0 


-v^  o 


“3  “  a  %2°s 


(45) 


IL 

-V. 
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Then  in  our  difference  approximation 

“3+  Kl  -  “3+  yu+  In  “4+  -  ^4+  In  “4+ 


V  ;  73+ \  + 


“3+  hi  "  “3+  h 


3n  t'sh- 


(A£) 


where 


(47) 


0o  =  0(x  =  0) 

0+1  =  0(x  +  6^). 

Expression  (46)  is  to  be  substituted  into  eq.  (43) . 

Note  that  J  is  not  assumed  continuous  across  the  boundary 
-  J_  -»•  0  for  6^  '♦■0,  6_  -*•  0,  but  for  a  finite  mesh,  we  do  not  use  the 
continuity  of  J  as  an  explicit  condition. 


D.  The  Constants  C. 


mn 


We  had 


(-1  (-1  _ P 

Cj^(x,u)»  2Tn  M(ulx,u)(i|i  \  L^Cv)  ^  (x,u)"'  \  [1™+ Ijl-y^cos^  "d? . 

^-1  -1  ^  Jo  (18) 


We  assume  that  within  a  given  medium  is  independent  of  x  and 
further,  that  M(m.1x,u)  is  Independent  of  u. 

Let 

lJ=  iiv  +  \ll-M.^  \ll-v^  cos  0 
We  can  expand 


(48) 


(49) 
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Then 


£ 


P^(V^)d0  - 


/■O 


P«(ii)P*(v)+2  2  P?(v)c08  in0  I  d0  - 
"  "  maO  “ 


•] 


n 


(50) 


The  second  eqviallty  arises  from  the  addition  theorem  for  Legendre 
Polynomials. 

We  use  also  the  Inverse  of  (47) : 


P^(n) 

and  define 

<|x*S 


(51) 


(52) 


The  differential  cross  section  as  a  function  of  the  cosine  cu 
of  the  scattering  angle  In  the  center-of-mass  system  is  expanded  as 


O  (u)  00 

«ys(u,aj)  - — 2  (2i+l)  f^(u)P^(a>), 

^“0 


(53) 


where  a-(u)  Is  the  total  macroscopic  scattering  cross  section  at  u, 

B 

from  which  It  follows  that  fQ(u)  ■■  1.  The  relation  between  v  and  cu 
Is  such  that 

ag(u,v)dv  -  ag(u,a))da.  (54) 
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Thus  we  can  write 


mn 


L^Cv)  (2i+l)f  (u)P.(a))da)2ir£  f  I  a  a  y? 

^X=0  ^  ^  t«0  p-0 


=  Y*  (2i+l)f  .(u)  2  b  2  a  2  ^  k<nSj. 
^i=0  ^  r=0  "^p=0 


(55) 


where 


C  l” 
‘^-1 


“^ipm  =  \  L  (v)v‘'P^((D)da>. 


For  hydrogen  the  lower  limit  of  the  integral  over  a>  is  zero. 
We  can  write 

^mO  “  T 


''m  =  i  L  <2i+l)fxJ 


i=0 


ilm 


(56) 


(57) 

(58) 


C„,  -  +  1  z  (21+l)f, 

2  i=0  *  ^  '£=0  '  ^ 


(59) 


The  sum  over  i  actually  will  be  cut  off  at  some  value  £  •*  L.  For  an 

2  1 

isotropic  flux,  <M,  >  =  and 


-m2  “  J  ^mO 


(59a) 


'^ipm  evaluated  recursively  in  the  code 


14 


E.  The  Calculation  Method 

We  use  three-point  difference  approximation  formulas;  for 
both  first  and  second  derivatives.  For  the  first  derivative  we  have 
a  weighted  average  of  the  forward  and  backward  derivatives,  thus: 


.  ( i  +  i  )  [i 
A  J 


^k+l"^k  .  1 
A 


A++  A_ 


A+  A_  \  A^ 

A“  "  a;;  j^k  ■  A“  "  ^-11  " 


(60) 


where 


M  |j^  “  M  ''k- 

^k  “  ^^“k^ 

\  “  \+l  “  \ 

A.  =  -  ^.i 

For  the  second  derivative  we  take  a  difference  quotient  between 
the  forward  and  backward  first  differences: 

i0 


F'^k+l'^k 

4-4-11 

.A  ^  J 

L  ^ 

A_  J 

’  a;?x[^  '*k+i  \  '“k-i  i 


(61) 


Suppose  is  the  flux  at  position  mesh  point  j  and  lethargy 
mesh  point  k.  Let  us  define 


15 


\  “  ■a^(a^+a_) 

2 

“  a^Ta++a_) 

X  -  2 

+  "  6+ <6^+  6J 


p,  s» 


■  A_(A^+  A_) 
2 

A_(A^+  A_) 

2 

B_(6^+  6 J 


(62) 


where 


*d+l  ■ 


®  x«  “  * 

«j  j 


Then 


M 

jk  “ 

\ 

^j,k+l  " 

(V 

'  ^jk  + 

"-'•j.k-i 

(63a) 

s 

jk 

^+1" 

+  fJt.)  0jic  +  Vi. 

•  ^j,k-l 

(63b) 

& 

jk 

^j+l,k  " 

i\+  K'^ 

'  <*jk* 

1 

1 

M 

>• 

K 

(63c) 

M 

2 

Ft 

^l+l,k' 

fit 

^ik"  ^i-l.k 

^  OK 

jk  ‘ 

*  8^+  6_ 

l7k+ 

5+ 

T'k- 

5_ 

”  ^*^16+  ^j+l,k  ‘  ■yk+'‘‘  ^-'>^k-^  ^^jk  ■*■  ^-'^k-^j-lA 

(63d) 


and  are  respectively  the  values  of  y  at  lethargy  k  in  the 
rlght-and  left-hand  Intervals. 
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The  Tj's,  ^i,»s  and  X's  are  functions  of  j  and  k.  Chough  we  omit 


the  indices  for  convenience.  Our  equations  become 


p=-l  qLl  ^P‘1  ~  ' 


(64) 


where  is  a  source  term,  if  any.  Here  the  SPM  differential  equation 
can  be  written 

_a°0  -  In  ■ 


where 


a 

In  H 

(65) 

0 

a 

=  a  -  Cs 

(66a) 

“  ^10 '^s 

(66b) 

a2 

“  7  ^2Q^s 

(66c) 

a3 

=  0  ^OZ'^s 

(66d) 

=  ^12 ‘^s 

(66e) 

7' 

0 

> 

1 

0 

11 

(66f) 

7^ 

1 

s 

a 

(66g) 

7^ 

=-'^t 

n 

(66h) 
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The  a’s  and  7's  are  constant  within  each  region'.  The  a  that 

pq 

appear  in  eq.  (71)  are  determined  as  follows: 

®+l.l  *  ■'‘+Vki^''k+  “k+l,+  ■'k+l,t  4tl,+> 


*0+1 


“k+1  “k+1  "  ^-1,+1  "®1,+1 


=  X. 


*+1,0  ”  "+  'yk+ 


4+  +  (\  + 1->  “k+  +  yl  4+5 


a 


0,0 


-c^  +  (\  +  n.)4  -  K  +  “k  -  M,o  '  ^1,( 


(67) 


If  we  are  not  at  a  boundary  point  and  if  the  spatial  integration 
mesh  is  constant  in  each  region,  then 

n,q  -  “-l,q 

To  solve  eq.  (64)  we  use  a  line  relaxation  method.  The  solution 
goes  as  follows:  Suppose  we  want  0jjj»  0  ^  j  ^  J,  0  ^  k  ^  K.  An 
initial  guess  ^j^^is  made  for  at  all  interior  mesh  points,  i.e., 

for  1  ^  j  ^  J"li  l^kj^K  -  1.  The  conditions  that  the  flux 
vanish  at  the  boundary  are  given  by  the  expressions 


"jK 


0, 


ok 


"Jk 


0  1  j  1  J 

1  1  k  ^  K 
1  ^  k  ^  k. 


(69a) 

(69b) 

(69c) 


0jQ  is  prescribed. 
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Given  an  estimate  of  the  flux  for  lethargies  (k-1)  and  (k+1) , 
we  can  get  a  new  estimate  for  lethargy  kby  the  relation  (64). 

Then  in  the  line  relaxation  procedure  we  get  from  eq.  (64) 


p=-l 


^  ^po  ^j+p,k  ^jk  ' 


(70) 


where 


^jk  -  Sjk 


1 

2 

p=-l 


®p,-l  ^j+p,k-l 


1 

-  2 


p=-l 


,jk  5j(i-l) 

P^i  ^j+P>k+l 


(71) 


We  have  put  indices  j,k  on  the  a-coef ficients  for  explicitness.  A 
superscript  (i)  represents  the  i—  estimate  of  the  flux.  We  start 
with  k=l.  ^j+p  k-1  then  just  the  prescribed  flux  at  u  =  0.  For 
k=2,  0^?^  0  is  the  initial  estimate.  Thus  g^Pcan  be  determined  from 

j  +p , z  j i 


eq.  (71) .  Solving  eq.  (70)  then  gives  0: 


(i) 


jl  • 


The  procedure  is  then 
s(i-), 


repeated  with  k=2,  etc.,  down  to  k  =  K-1.  We  now  have  0jj^' everywhere . 
The  process  can  clearly  be  iterated  to  give  successively  higher 
approximations.  It  can  be  shown  that  if  the  process  converges,  it 
converges  to  the  correct  answer. 

The  relaxation  procedure  could  as  easily  have  been  carried  out 
along  lines  of  constant  x.  The  equations  that  then  replace  (70)  and 
(71)  are 

=  f(i) 


1 

2 

q=-l 


aJk  jj(i)  ^  fv 

®oq  ^j,k+q  ^jk 


f(i)_  s  -  7  a^^  -  2  aJ^ 

-  ^jk  ^l,q  ‘^j+l,k+q  J_^^-l,q<^j-l,k+q 


(72) 

(73) 
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We  use  row  Iteration,  although  a  column  iteration  or  alternate 
row  and  column  iteration  are  other  possible  procedures.  Eqs.  (70) 
and  (72)  are  of  the  form, 


Vn-l  +  '=n  +  'n  K-tl  '  %  ,  0  <  n  <  N. 


(74) 


Eq.  (74)  is  solved  by  use  of  two  auxiliary  parameters  and 
which  are  defined  recursively  as 

c. 


n 

b_+  a^Z 


n 


n^'n-l 


^n~^n  ^n-1 
^n'^  ^n  ^n-1 


(75) 

(76) 


with 


0. 


(77) 


Then 

<'n  -  fn  +  'o+l' 

Thus  we  have  P  ,  Z„.  Given  P_  i ,  Z_  i ,  we  can  compute  P_,  Z  by 
o’  o  n-i*  n-l’  n’  n  ■' 

eqs.  (75)  and  (76)  out  to  n  =  N-l.  is  known.  Given  we 

can  compute  0^  by  eq.  (78),  back  to  n*»l.  Thus  the  complete  solution 
is  given. 

In  a  spatially  symmetric  situation,  we  take  J  0  along  the 

plane  of  symmetry,  with  -J  <  j  ^  J.  Symmetry  implies  that  0-^  0_j^  ^ 

and  that  a.  =  a  i  .  For  relaxation  along  lines  of  constant 
i,q  “•Ljtl 

lethargy,  we  have  to  solve  equations  of  the  form  (74)  for  0  <  n  ^  N 
with  the  added  relation 

^o’^-l  +  Vo  +  '=o’^l  “ 
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But  because  of  the  symmetry,  and  Thus 

2a. 


Qo 


(80) 


Comparison  with  eq.  (78)  gives 


o 


= 


(81) 

(82) 


can  then  be  found  from  eqs.  (75),  (76),  and  (78)  for  0  ^  n  <  N. 
Relaxation  along  lines  of  constant  x  can  be  carried  out  either 
from  N  to  0  from  0  to  N.  However,  for  j  =  0,  eq.  (23)  becomes 


^ok  “  ?ok 


h,k+q‘ 


..th 


The  i—  estimate  for  0^^  is  used  in  proceeding  from  K  to  0. 


(83) 


The 


s  t» 

(i-1) —  estimate  is  perforce  used  if  we  proceed  from  0  to  K. 
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F.  Source  Conditions 

Assume  we  have  a  plane  isotropic  source  5(x)5(u)  at  u  =  0. 

For  a  uniform  medium  the-  slowing  down  density  at  u  =  0  can  then  be 
written 

q(x)  =  Y  (a|xj),  (84) 

since  everything  that  scatters  slows  down  at  the  collision  point. 
a  is  the  cross  section  at  u  =  0.  If  there  is  little  absorption  near 
the  source  energy,  q(x)  is  nearly  constant  except  near  x  »  0.  But 
except  immediately  around  u  =  0,  we  have 

q  =  (85) 

where  ^  is  the  mean  lethargy  increase  per  collision. 

Thus  the  flux  is  given  correctly  several  ^-units  and  more  from  the 
boundary  if  we  take 

0(x,O)  [xj  ).  (86) 

We  use  this  as  a  boundary  condition.  It  must  be  stressed  that  0(x,O) 
is  not  a  valid  approximation  to  the  flux  very  close  to  u  =  0.  It  is 
the  condition  that  must  be  used  with  the  SPM  equation  to  give  the  flux 
correctly  away  from  u  =  0.  With  this  normalization  both  SPM  and  age 
theory  reduce  to  the  same  asymptotic  result  for  large  u,  as  they 
should. 

At  an  interface  x^  we  get  the  condition  by  assuming  that  the 
number  of  neutrons  slowed  down  per  unit  time  from  u=o  between  x^-6B_ 
and  is  given  correctly  by  assuming  that  the  flux  at  u  »  0 
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is  continuous  and  linear  in  the  intervals  (x  -e6  ,x  )  and(x  ,x  +  eo). 

6^'  and  6_  are  the  distances  from  to  the  next  mesh  points  to  the 
right  and  left,  respectively.  The  natural  requirement  would  be  to  let 
e  =  1,  but  in  order  to  get  an  expression  which  reduces  to  the  correct 
one  when  the  media  on  the  two  sides  of  the  interface  are  identical, 
we  let  e  approach  zero.  The  precise  condition  at  a  single  point  is 
not  especially  critical,  in  any  case.  Thus  we  write  for  the  slowing 
down  density 

o,  . 

<l(x)  =  -j-  Ej^Cjs  +  (x-XQ)j)  Xq  <  X  <  x^  +  e6^ 

a_ 

=  Ei(|s  -  a_(xQ-x)  |)  ,  <  x  <  x^  ,  (87) 

where 

0^  =  total  macroscopic  cross  section  in  interval  x  <  x  <  x^+€5_j_ 

a_  =  total  macroscopic  cross  section  in  interval  XQ-e6_<  x  <  x^ 

s  =  total  number  of  mean  free  paths  between  the  origin  and  x^ 

We  assume  also  that 

q(x)  =  ,  Xq  <  X  <  x^  +  €6^ 

= 

Here 

=  mean  lethargy  increase  per  collision  in  interval  X^<X<Xq+€5^ 
=  mean  lethargy  increase  per  collision  in  interval  Xq-€5_<x<x^ 
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Then 


\  q (x)  dx  -  \  7“  ^1  (  h"^^-  I  ^  ‘^^■*’  \  ^1  ^  1®'*''^+  i  ^ 


rS+€0+6+ 

7  \  Ej^(lylf)dy. 

'^s-ea  6 


Also 

J  X^-€6 


q(x)dx  *»  ^_a 


^  X, 


0(x)dx  +  ^ 


•€6 


po+e6+ 

+‘^+  y  ' 

V 


0(x)dx. 


(90) 


Now  in  the  limit  €  ->  0,  we  get  by  equating  expressions  (89)  and  (90) 


Y  £3^(5)  (€0^6^+  €o_6j  =  0(Xq)  (€|^0^6^  +  e^_a.6.)  , 

or  finally 

0(Xo)  «  7  Ej^(s)  4_a_5_  * 

For  an  interior  point,  this  reduces  to 


(91) 


0(Xo)  -  Jej^(s),  (92) 

as  expected. 

The  source  gives  a  problem  in  a  numerical  calculation 
because  it  is  singular  at  the  origin.  The  criterion  we  use  to  avoid 
the  singularity  is  that  the  total  number  of  neutrons  slowed  down  per 
unit  time  from  u  ■  0  in  the  interval  -6_<  x  <  6^  is  given  by  a 
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three-point  formula  in  terms  of  the  continuous  and  piecewise  linear 
flux  Sj(x)  at  u  =  0. 

Thus  eqs.  (87)  and  (88)  hold,  with  =  0,  s  =  0,  e  =  1. 
f' ^  a  o 

1  q(x)dx  ={  El  (  ia_x|)dx +\  -5^  E,  (  o.x  )dx 

J-5.  J-5.  Jo 


r  c^.5. 

\  ,lV  1-E2(a,.; 

■j  Ej^(y)dy+Y\  Ej_(y)dy  «  - ^ - 

Cl  w/  0 


1 

7 


■a  6 


°+^+ 


e  +  cT_6_E2^(a_6_) -e  +  o^S^Ej^  (a^6^)J 


We  have  used  the  relations 
*00 


I 


(93) 


E]LCy)<iy  =  E2(p) 

(94) 

£2(0)  =  1 

(95) 

-P 

(96) 

E2(p)  =  e  -p  Ej^(p) 

Also, 
i5 


^  q(x)dx 


o 

a_0(x)  dx  +  \  5^a^0(x)  dx 

6  o 


=  ^  l0(-6.)  +  0(0)]  +  40^  ^j0(O)+  0(6^^ 


(97) 
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Since  by  assumption,  -B_  and  5^  are  not  interface  points. 


and  so  equating  the  right-hand  sides  of  eqs.  (93)  and  (97)  and 
using  (98) ,  we  have 


•o  6 


-‘^+s+ 


Y  [2  -  e  +  a_6_E^(a_6_)-e  +  o^6^Ej^(0^8^^  = 

=  5^.8^)  0(0)  +  I  a_6.Ej^(a.6_)  +  J  a+8^E^(a^6^) 


or  finally. 


-0  6 


~°+\  1 


0(0)  = 


2  -  e  +  ^  8_6.Ej^(a_6_)-a  +7  (0^6^) 

+  40^6^ 

Eq,  (99)  gives  the  value  we  use  for  u  =  0  at  the  origin. 


(99) 
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G.  Boundary  Conditions 

Since  we  are  working  with  an  equation: essentially  of  the 

\ 

elliptic  type,  we  specify  one  boundary  condition  all  around  the 
boundary.  The  computations  use  a  rectangular  region  (a<x<b,  0<u<u^; . 

The  source  condition  at  u=0  has  been  discussed.  V7e  use  the  condition 
0(x, u^)  =  0  at  the  cutoff  lethargy  u^.  It  was  shown  in  the  previous 
summary  report^  that  the  precise  condition  at  u^  does  not  effect 
the  solution  except  in  the  immediate  vicinity  of  u^,  if  tne  condition 
is  at  all  a  reasonable  one.  The  present  spatial  boundary  conditions 
are  that  the  flux  vanishes  at  a  and  b.  One  could  use  instead 
derivative  boundary  conditions,  i.e.,  that  is  equal  to  some 
constant  on  the  boundaries.  This  is,  in  fact,  the  condition  on  the  line 
of  symmetry  for  the  symmetrical  case,  where  '^  =  0. 

H.  Inelastic  Scattering 

Inelastic  scattering  must  be  treated  as  an  absorption  followed 
by  a  reemission  at  a  lower  energy,  as  in  multigroup  theory.  It 
presents  two  problems,  one  intrinsic  and  one  purely  practical.  The 
practical  problem  is  that  of  time  and  memory  requirements  in  the 
machine  computation.  We  get  an  effective  source  of  strength 

S(x,u)  =  ^  a(x;  u',  u)0(x,u')du'  ,  (100) 

where 

a(x;u',u)  “  inelastic  macroscopic  scattering  cross  section  at  x 
for  neutron  of  lethargy  u'  to  scatter  to  u,  per  unit 
lethargy  range  about  u. 
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S  is  treated  as  an  inhomogeneous  neutron  source.  The  appropriate 
term  in  the  SPM  equation  is  the  first  collision  density  of  neutrons 
Whose  source  strength  is  given  by  S,  , 


l('o(x,u)  = 


1 

7 


where 


Ej 


a(u)  \  E,(|s|)  S(x',u)dx' 


(101). 


s  =  y  a(u,x")dx'', 

\J  x' 


Just  as  in  multigroup  diffusion  calculations,  however,  we  will  not 
take  account  of  the  transport  of  these  neutrons  between  the  point  of 
inelastic  scattering  and  the  point  at  which  they  next  collide.  That  is, 
we  take 


=  S(x,u).  (102) 

This  approximation  eases  the  very  considerable  computation  that  would 
be  required  to  use  eq.  (101) . 

The  intrinsic  problem  with  inelastic  scattering  lies  in  the  fact 
that  any  source  in  the  interior  of  the  u-interval,  i.e.,  for  u  >  0, 
can  propagate  in  both  directions  in  lethargy.  While  the  propagation 
to  lower  lethargies  is  damped  fairly  rapidly  with  decreasing  u,  this 
part  of  the  contribution  is  certainly  unphysical.  Further,  the  flux 

i 

is  not  correct  even  above  this  interior  source  lethargy  until  about 
an  interval  Au  5  away.  It  is  certainly  true,  however,  that  if  the 
flux  inelastically  scattered  to  an  interval  whose  width  is  of  the 
order  of  ^  is  small  compared  to  that  slowed  dovn  by  elastic  scattering, 
the  error  is  small.  That  is,  in  that  interval  the  error  is  small  by 
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hypothesis,  and  in  any  case  does  not  propagate  very  far.  The  latter 
property  implies  that  the  error  is  not  cumulative,  but  only  local 
so  that  it  is  limited  in  magnitude  everywhere. 

It  should  be  remarked  that  these  difficulties  are  present  in 
part  in  ordinary  multigroup  theory  whenever  the  relation. q  =  ^o0  is 
used  to  connect  the  slowing  down  density  and  the  flux. 

The  code  takes  inelastic  cross  sections  in  two  forms:  either 
as  the  cross  section  a£^(u»,E)  at  u*  per  unit  final  energy  about  E, 
tabulated  as  a  function  of  u'  and  E,  or  as  an  excitation  cross  section 
f^j^CuSQ^)  tabulated  as  a  function  of  u'  and  the  excitation  energy  Q^. 

In  the  available  tabulated  inelastic  cross  section  data 
be  available  as  interpolated  in  the  initial  lethargy  directly  from 
experiment,  for  the  lowest  excitation  levels  Q^.  At  higher  incident 
energies  one  resorts  to  calculated  cross  sections  in  the  fom  a^^(u',E). 

In  general,  precisely  to  the  total  observed 

inelastic  cross  section  at  u'  because  of  experimental  inaccuracies. 

We  therefore  define  the  normalized  excitation  cross  section 


gi(u)  =  fin(u,Qi) 


(103) 


gj^  has  the  required  property  that 

I  gi(u)  =  •  (10^> 

a^^(u)  is  the  cross  section  for  neutron  production  by  inelastic 
scattering.  In  terms  of  the  (n,  n')  and  (2,  2n)  cross  sections  it  is 


-  “n,n. 


(105) 
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Since  .E)  is  computed  while  is  measured, 

integration  of  a^^(u,E)  over  E  will  not  in  general  give  tT^j^(u). 
Again  We  must  impose  a  normalization.  In  the  code,  it  is  asstomed 
that  0.  (u' ,E)  is  linear  in  E  between  each  successive  tabulated 
values.  It  is  further  assumed  that  if  E  takes  on  the  values  Ej|^ 
and  if  there  are  F  values  of  i  in  the  tabulation^  that 


®P+1  “  ° 

Oin(u'>0)  * 

^i+1  ^i 


0 


^in(u' «E)  =0,  E  >  E^. 

Then  we  can  get  a  normalized  cross  section 

o,-_(u’) 

Oin'(u',E)  =  a^^(u',E)  — ^ -  (106) 

t 

In  terms  of  gj^  and  the  differential  inelastic  neutron 

cross  section  is 

a(u',u)  =  E  cT^^(u' ,E)  ,  u'<  u^ 

* 

=  2  g.(u')  6  (E'-E-Q.)E,  u'  ^  u,  (107) 

i  ^  ^  h 


Here  u, 
*^1 


is  the  smallest  lethargy  at  which  given. 


For 


Uj^^  <  u'  <  »  where  Uj^  is  the  largest  lethargy  for  which  we  use 

the  (u'  ,E)  data,  the  codes  takes 


'tn 


(108) 


Eq.  (107)  asstimes  that  the  excitation  energy  Qj.  equals  the  neutron 

energy  loss  in  the  c.m.  system.  While  not  strictly  true,  this  is 
an  adequate  approximation  for  our  purposes. 


30 


I 

■'-4 


Elsewhere  we  interpolate  linearly  in  u*,  for  both  a£^(u',E)  and 

gi(u’). 

In  eq.  (107), 


Note  that 


E«  =  E^e 


-u 


(1C9; 


Qi  <  S^e 

Substituting  eq.  (107)  in  eq.  (100)  and  using  eq.  (107),  we 

have,  dropping  the  spatial  variable  x, 

>  u 


^  1  "  r ' 

^  i.  y 


S(u)  = 


,u)0(u')du' 


=  (.(U. 

w'  o 

=  r  .(u., 

■J  O 

=  \  (u',E)E  0(u')du'  +  E  E  \  g^(u')  6  (E» -E-Q^)  0(u')  du,' 

^  '^u,. 


E)E  0(u»)du', 


U  >  U,  (111) 
‘^1 

In  our  mesh  scheme,  the  lethargy  takes  on  the  value  Uj^,  and 
integration  is  by  the  trapezoidal  rule.  Let  us  define 

Gk  =  -2^  V'*'  ^(^i+P‘^in^“i+l'^k^}  ^^i+l'^i^  ' 


with 


i=l 


ko  -  k  . 


(112) 


“k  <  “k. 


=  k,  , 


''k  ^  ^k. 


(113) 
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'^k  related  by  the  standard  expression 


E 


k 


-u 


k 


Then 


i: 


,E)E0(u  )du  =  Gi,  . 


Also, 


r'\ 

)  Si(u': 

^  II 


E  \  g^(u  )6(E'-Ej^-Q^)0(u')du'  = 
u 


g^(in  ^  )  0(in  ^  )6(E'-E, -Q.)  ^ 
^  E  E  ^  ^  E 


=  8^(in  )  0(in 

H  is  again  the  step  function  defined  previously: 


H(e)  =  1, 

=  0, 


e  >  0 

i  <  0. 


Let 


and 


a. 


ik 


Ek+Qi 


=  1  +  yr-  e 


Vik  =  in  a, 


ik' 


(i 


(115) 


(3) 

(116) 


(117) 
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Then 


(118) 


Substituting  eqs.(114)  and  (115)  in  eq. (Ill)  and  using  eqs. 

(116)  and  (118) ,  we  get  finally 

S(uj^)  =  +  ^  8i(V''ik)^^<V'^ik)«(\-W- 

Eqs.  (113) , (109a) ,  (3),  (116),  and  (117)  define  all  the  parameters. 

Eq.  (119)  holds  at  every  spatial  mesh  point j  in  the  inelaacic  scatteringy 
medium. 

To  determine  what  value  of  S(uj^)  to  use  at  a  boundary  point  of 
an  inelastic  scattering  medium,  let  us  consider  a  situation  where 
there  is  one  inelastic  scattering  medium  to  the  right  of  spatial 
mesh  point  j  and  another  to  the  left.  Then  we  assume  that  the  total 
source  of  inelastically  scattered  neutrons  is  represented  as  being, 
continuous  at  Xj  and  linear  between  and  and  between  Xj  and 


'j+1* 


Thus 


r  J+i  6_r 

\  ,  S(x,Uj^)dx  =  -z-LSj-i^k  +  S. 

jk]'^  3r  hk"^  ®j+Uk]  “  ^^''k 

(120) 

Here 

Sjk  =  S(Xj,Uj^)  . 

(121) 

and 

Sj+l,k“  ^^“kp* 

(122) 
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Then 


6,S(u,  ,)  +  6  S(u,  ) 

T  iC-r'  -  '  Ic- 


6_^  +  6 


(123) 


In  the  spirit  of  the  derivation  of  the  flux  at  an  interface  for  u-----.';. 
eq. (123)  should  hold  not  for  the  density  S  of  inclastically  scattere. 
neutrons,  but  for  the  cross  sections  0^^(u,E)  and  g^(u).  However, 
prescription  (123)  is  much  simpler,  in  addition,  the  overall  error 
is  small,  probably  much  smaller  than  the  error  involved  in  the 
assumption  (102) . 
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III. 


Results 


A.  Carbon 

Comparisons  were  made  of  SPM  results  with  those  of  age 

2 

theory  and  of  a  method  due  to  Kaper.  The  medium  was  carbon  with 
an  assumed  constant  cross  section.  There  Is  no  absorption. 

The  flux  given  by  age  theory  can  be  written 

0(x,  u)  = 

Here  the  geometry  Is  that  of  an  Infinite  plane  source  of  strength 
S(x,  u)  =  6(x)  6(u) 
with  no  first  flight  correction. 

T  Is  the  neutron  age,  given  by 


The  age  theory  solution  Is  tabulated  In  Table  I. 

2 

Kaper  has  developed  a  treatment  which  Is  essentially  a 
consistent  expansion  In  powers  of  1/A,  valid  also  for  variable 
cross  sections.  Since  one  can  look  on  the  SPM  as  an  expansion  In 
powers  of  1/A,  though  not  a  consistent  one  (since  It  gives  certain 
higher  order  terms) ,  It  Is  of  Interest  to  compare  the  two  methods. 
Kaper  has  given  results  for  the  first  order  solution;  l.e. ,  that 
good  to  linear  terms  In  1/A.  The  SPM  equation  contains  In  addition 
the  terms  In  1/A  and  certain,  but  not  all,  higher  terms.  For 
carbon  (A  =  12) ,  the  differences  should  not  be  great. 
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For  the  comparison  we  use  Kaper's  solution  including 
first  flight,  that  is,  with  a  first  flight  correction,  normalized 
to  the  age  theory  results  at  Zx  =  1,  u  =  15.  Of  the  points  given 
by  Kaper,  this  is  the  one  for  which  age  theory  should  be  best 
(small  X,  large  u) ,  and  so  the  results  were  normalized  accordingly. 
This  solution  is  shown  in  Table  II. 

We  have  two  sets  of  SPM  results  for  this  problem.  Both 
used  an  Ej^  source.  One,  given  in  Table  III,  is  the  solution  of 
the  third-order  equation.  The  other,  shown  in  Table  IV,  is  the 
solution  when  the  third-order  terms  (i.e.,  7^  and  7^)  are  assumed 
zero.  Note  that  the  effect  of  the  third- order  terms  in  this 
problem  is  negligible. 

In  comparing  the  solutions  given  in  Tables  I-III,  we 
note  that  agreement  is  very  close  at  smaller  distances,  but  for 
larger  values  of  x,  the  SPM  solution  is  flatter  -  smaller  for  large  u 
and  much  larger  for  small  u  -  than  either  the  age  or  Kaper's  solu¬ 
tion.  Certainly  the  Increase  for  small  u  is  an  improvement.  The 
differences  for  large  u  are  not  great  and  it  is  not  clear  whether 
or  not  they  are  significant. 
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.TABLE  I.  DIFFUSION  THEORY  (A  -  12) 


-.1 

®  7  \|tu  ® 


1 

2 

5 

8 

10 

1 

1.06 

0.726 

0.053 

- 

- 

3 

0.663 

0.585 

0.244 

0.048 

0.010 

5 

0.522 

0.484 

0.286 

0.108 

0.044 

8 

0.417 

0.398 

0.286 

0.156 

0.089 

10 

0.374 

0.360 

0.277 

0.170 

0.108 

12 

0.342 

0.331 

0.266 

0.177 

0.122 

15 

0.306 

0.299 

0.251 

0.181 

0.134 
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TABLE  III.  SPM  FLUX  (A  -  12)  -  WITH  THIRD  ORDER  TERMS 


1 

2 

5 

8 

10 

1 

1.07 

0.760 

0.098 

- 

- 

3 

0.676 

0.605 

0.280 

0.070 

0.020 

5 

0.534 

0.499 

0.313 

0. 131 

0.053 

8 

0.426 

0.408 

0.302 

0. 165 

0.079 

10 

0.381 

0.368 

0.286 

0.167 

0.083 

12 

0.347 

0.336 

0.267 

0. 162 

0.082 

15 

0.305 

0.296 

0.240 

0.149 

0.076 
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TABLE  IV.  SPM  FLUX  (A  -  12)  -  WITHOUT  THIRD  ORDER  TERMS 


1 

2 

5 

8 

10 

1 

1.09 

0.728 

0.096 

- 

- 

3 

0.691 

0.611 

0. 270 

0.069 

0.020 

5 

0.541 

0.504 

0.309 

0.127 

0.051 

8 

0.430 

0.301 

0.162 

0.077 

0.384 

0.370 

0.285 

0.165 

0.082 

12 

0.347 

0.337 

0.267 

0.161 

0.081 

15 

0.306 

1 

0.297 

0.241 

0.149 

0.076 
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B.  Hydrogen 

A  nuiid}er  of  problems  were  r\m  for  hydrogen,  since  the 
ability  to  treat  hydrogen  was  considered  a  major  part  of  the 
program.  Some  used  an  IS-Mev  source  and  some  an  8-Mev  sotirce, 
with  X-  and  u-meshes  of  various  sizes.  Results  for  8-Mev  sources 
are  not  yet  available.  All  of  the  18-Mev  problems  failed.  The 
reason  Is  important  and  instructive. 

All  of  the  18-Mev  problems  Included  the  third  order 
terms.  Basically,  they  failed  because  and  a_j^  (eq.  67) 

were  too  large.  a^Q  could  therefore  be  small  and  even  negative. 

If  we  were  dealing  with  a  delta  function  sotirce,  for  Instance, 
i.e. ,  a  source  of  unity  at  (j,k)  ■'  (0,0),  where  j  and  k  are 
respectively  the  x  and  u- indices,  this  would  imply  that  0(1,O)<O. 
This  is  obviously  inadmissible.  The  situation  with  the  Ej,  source 
was  somewhat  more  complicated,  but  fundamentally  the  same.  Not 
only  were  some  fluxes  negative,  but  the  error  propagated  from 
iteration  to  iteration  and  the  iteration  procedure  rapidly 
diverged.  It  should  be  remarked  that  a^^  and  a^^  can  become 
too  large  only  in  the  case  of  variable  cross  sections. 

The  conclusion  from  all  this  is  that  if  the  cross  section 
variation  is  too  great,  trouble  can  be  expected  from  the  third- 
order  terms.  This  bears  on  the  problem  of  treating  resonances. 

If  we  had  hoped  to  Improve  on  the  second  order  equation  by  adding 
third-order  terms,  one  must  be  extremely  careful  with  variable 
cross  sections.  We  have  not  been  able  to  deduce  any  real  rules 
for  predicting  when  the  third-order  terms  help. 
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C.  Modification  of  Diffusion  theory  Current  by  SPM 

To  see  how  the  additional  terms  in  the  SPM  equation 
modify  the  diffusion  theory  expression  for  the  current,  we  have 
two  procedures.  In  one,  we  can  assume  that  the  change  from  age- 
diffusion  theory  is  a  small  perturbation  and  apply  the  derivatives 
in  the  additional  terms  to  the  age-theory  solution.  In  the  other, 
we  look  at  the  ratio  of  the  additional  terms  in  the  difference 
expression  for  the  ciarrent  to  the  diffusion  term  for  one  of  our 
numerical  solutions.  We  will  consider  both  methods,  for  a  medium 
with  constant  cross  section  and  no  absorption. 

The  age-diffusion  solution  is 


with 


u 


The  SPM  expression  for  the  ctirrent  in  a  constant  cross 
section  medium  is 


J  “  ■'S  [•‘sH  -  +  T'5°'3> 

We  are  then  interested  in  R,  the  ratio  of  the  additional  term 
in  this  expression  to  the  diffusion  term. 


R 


®<3 
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In  our  case 


wLTsi .  n  +  ^ .  s  + 

"<3  ^02  ^ 

for  an  Isotropic  flux. 

We  take  derivatives  of  the  unperturbed  age  theory 

solution 


in  0  “  in  (oO  -  y  in  (Att)  -  y  in  t  -  ^ 
F  ^  7F 


Thus  the  analytical  approximation  to  R  Is 


For  carbon,  In  our  problem, 

4  =  0.15779 

=  75  "  0.0A953 

I-Vq  o 


Then 


R 


an 


0.10826 


0.078895 

■3r 


x2] 


A3 


We  compare  this  with  the  nvimerical  approximation, 
in  which  the  derivatives  are  evaluated  by  differences  from  our 
numerical  results.  The  results  are  shown  in  Table  V. 

TABLE  V.  EFFECT  OF  HIGHER,  ORDER  TERMS  ON  THE  CURRENT 


X  (mf  p) 

u 

^anal 

D 

num 

2 

1 

0.111 

0.108 

2 

0.0683 

0.0667 

2 

5 

0. 0304 

0.0299 

2 

10 

0.0157 

0.0098 

5 

5 

0.0197 

0.0194 

5 

0.0130 

0.0075 

10 

5 

-0.0188 

-0.0201 

10 

10 

0. 0034 

0.0022 

We  note  that  except  near  the  origin  and  for  small 
lethargies  the  SPM  correction  to  the  current  is  small  -  of  the 
order  of  a  few  percent  or  less.  Where  and  differ 

appreciably,  the  correction  is  negligible. 
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IV.  Operating  Instructions  -  SPM  Code 

( 

— 

A.  Changing  Tape  Designations 

The  3  magnetic  tapes  required  by  the  SPM  code  are  assigned 
by  a  small  subroutine  called  CHANTP.  The  copy  of  the  SPM  deck 
Issued  at  the  time  of  this  report  defines  the  problem  Input  data 
tape  to  be  logical  tape  number  2,  the  output  tape  to  be  logical 
tape  nvonber  3,  and  the  tape  which  stores  complete  nuclide  data 
for  a  series  of  nuclides  to  be  logical  tape  number  5. 

The  tape  designations  may  be  inconsistent  with  the  monitor 
system  of  a  given  7090  Installation.  Therefore  a  simple  process 
for  making  the  necessary  changes  has  been  provided.  The  subroutine 
CHANTP  Is  listed  below. 


SUBROUTINE  CHANTP 

(Conmon,  dimension,  and  equivalence  statements) 
NUTAPE  =  5 
KTIN  =  2 
KTOUT  =  3 
RETURN 

END  (1,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0) 


Statements  of  CHANTP  subroutine 

KTIN  defines  the  problem  input  tape,  KTOUT  indicates  the 
output  tape,  and  NUTAPE  denotes  the  nuclide  data  tape.  To  redefine 
tapes  for  SPM: 

1.  Punch  the  proper  logical  tape  number  In  the  3  statements 
which  define  KTIN,  KTOUT,  &  NUTAPE, 

2.  Remove  the  old  binary  CHANTP  program  and 

3.  Recompile  CHANTP  (and  possibly  execute  a  SPM  problem). 

The  use  of  the  new  binary  CHANTP  subroutine  in  the  SPM  binary 
deck  or  tape  makes  the  code  usable  at  the  installation. 
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B.  Operation  Instructions  for  a  Problem  Run 

1.  Since  most  7090  installations  will  accept  Fortran  codes 
only  when  under  monitor  control,  only  one  provision  (see  below) 

for  console  control  of  the  SPM  code  has  been  made.  It  is  generally 
necessary  to  provide  the  7090  machine  operator  with  only 

a.  the  SPM  code  in  card  deck  or  magnetic  tape  form, 

b.  the  input  data,  usually  on  cards,  for  the  series 
of  problems  to  be  run, 

c.  the  nuclide  data  storage  tape  (NUTAPE)  which 
generally  will  not  be  written  on. 

Sense  Switch  6  must  be  up.  In  the  usual  case  nothing  further 
can  be  done  by  the  user  of  the  program. 

2.  Console  Control  of  the  SPM  code 

If  sense  switch  6  is  depressed  the  code  will  stop  after 
computing  and  outputting  a  complete  array  of  fluxes.  The  iteration 
number  of  the  flux  computation  loop  will  be  printed  both  on-and-off 
line.  This  console  operation  should  be  used  only  when  time  con¬ 
siderations  indicate  that  the  problem  being  run  should  be  terminated. 
Depressing  sense  switch  6  will  allow  the  program  user  to  obtain  all 
information  available  concerning  the  last  flux  iteration  at  the 
time  the  problem  was  stopped.  The  program  user  may  then  alter  the 
input  of  convergence  conditions  so  that  the  further  use  of  sense 
switch  6  will  become  unnecessary. 

NOTE:  In  the  following  input  table  we  use  the  following  definitions: 
^1  *  ®-l,-l  ^6  ®1,0 

^7  “  ®-l,l 
^8  “  %,1 

^4  “  ®-l,0  *^9  “  ®1,1 

S  ■  “o.O  -*10  -  Sjk 

where  the  a^^  and  are  those  in  eq.  (64). 


46 


“2  ■  ^0,-1 
‘'3  ‘ 


1 


The  columns  described  refer  to  the  72  columns  available  on  a 
standard  IBM  punched  card. 


Input  Data  for  SPM  Code 


Category  No. 

General  : 

Description 

Specific  Input  Quantity 

1. 

i 

i 

Titles  and  ;  (12A6) 

problem  ti 

‘  Identification 
;  Data 

j 

:  j 

_ _ i _ 

Date  in  columns  (1-12) 

A  problem  identification 
number  in  columns  (13-18) 

Any  comment  in  columns  (19-72) 
This  data  is  printed  in 
output  and  as  page  headings. 

^  1  ^  ^  I 

2.  j  Fixed  point  :  (2AI3)  j 

1  or  integer  i 

input  data: 

Note  that 

1  blank  columns  ! 

J  are  treated 

1  as  zeroes, 

j  Hence  all 

integer  data 

1  must  occupy  : 

1  right  hand 

1  columns  of 

their  fields 

MXIFL, 

Maximum  number  of  iterations 
before  problem  terminates 
itself.  It  is  placed  on  the 
i  first  card  of  category  2 
Columns  (1-13) 

i 

i  KLIMIT, 

I  Maximum  number  of  lethargy 

1  levels  included  in  convergence 
tests.  Columns  (4-6) 

j  J2, 


the  number  of  spatial  points. 
Including  boundaries. 

Columns  (7-9) 


the  number  of  lethargy  points. 
Including  boundaries. 

Columns  (10-12) 


IPCY, 

j  the  number  of  iterations  on 
the  flux  per  printout  of 
abbreviated  convergence  data 
(see  output  description) . 
Columns  (13-15) 

J 

I  IMALP, 

i  the  number  of  flux  iterations 
at  which  an  indication  of 
"Excessive  Iterations"  is 
printed  out  both  on-line  and 
off-line  and  beyond  which  all 
flux  convergence  data  is 


Category  No. 

General 

Description 

Format 

Specific  Input  Quantity 

2. 

printed  but;  the  problem 
continues,  however. 

1  Columns  (16-18) 

i 

(  ; 

1 

1 

MKO, 

the  mesh  size  of  the  lethargy 
output  printout.  Though  the 
computation  is  carried  out  in 

1  its  latest  stages  for  all  of 
i  the  K2  lethargy  points,  output 
;  occufs  for  only  those  points 
i  R=l,  1+MKO,  ...  1+nMKO,  ...  K2. 
Columns  (19-21) 

i  . 

1 

i 

[  ! 

1  1 

f  j 

i 

i 

j 

1 

MJO, 

the  mesh  size  of  spatial  point 
output  printout .  Only  the  fluxes 
corresponding  to  j=l,  1+MJO,  ... 
1+nMJO,  ...  J2  are  printed  out. 
Columns  (22-24) 

1 

i 

j 

f 

I 

[ 

( 

1 

j 

j, 

1 

IDEBUG, 
an  indicator . 

;  If  IDEBUG=20,  there  is  no  nuclide 

1  data  storage  tape.  Option  5, 
category  input  then  computes  or 
obtains  the  necessary  nuclide 
data. 

Columns  (25-27) 

1 

1 

|- 

1 

f 

,  MMT, 

the  number  of  entries  in  the  mesh- 
size  table  MIAB  (see  below) . 
Columns  (28-30) 
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MTAB(I),  1*1, MMT. 

A  table  of  the  mesh  sizes  to  be 
used  successive  in  the  compute- 
flux-array  computation.  The 
table  consists  of  pairs  of  mesh 
sizes  for  space  pt.  mesh>  and 
lethargy  pt.  mesh,  respectively. 
Coarse  mesh  pairs,  monotonically 
decreasing,  precede  the  finest 
mesh  size  pair  (l,l)whlch  are  then 
followed  by  two  zeroes.  By  fair 
the  most  usual  input  for  this 

Category  No. 


2. 


General 

Description 


p 

( 


I 

I 

i 


Format 


Specific  Input  Quantity 


;  table  of  MTAB,  where  MMT-A, 
I  1,  1,  0,  0.  NextCMMT  x  3 
I  columns) 


Is 


MIFLAG, 

the  number  of  entries  for  the 
:  Input  table  of  flux  arguments 
'  data  (see  table  ALPLAY  In 
category  3  Input) .  Next  three 
columns . 

!  I  KBGIHS, 

I  I  the  number  of  the  lethargy  point 

I  I  at  which  the  arbitrary  Input  of 

I  li  the  DIO  (Inhomogeneous  source 

j  I  terms)  begins, 

i  f  Next  3  columns. 

I  4' 


I  KIST, 

I  the  number  of  the  last  lethargy 
I  pt.  at  which  there  will  be 
I  arbitrary  Input  of  the  D-,q 
I  terms.  All  lethargy  points 
from  KBGIHS  to  KIST  will  have 
I  this  Input  whenever  KBGIHS  /  0. 
I  Next  3  columns. 


I  INONU, 

[  the  number  of  geometrically 
f  defined  regions,  each  of  con- 
I  stant  nuclear  composition. 

I  Next  3  columns. 


[  ID  10  FS, 

f  If  this  indicator  =  0,  fission 
Is  treated  as  absorption.  At 
I  this  time  the  indicator  must  be 
!  zero. 

I  Next  3  columns. 


I  ID  10  IE, 

i  If  this  indicator  »  0,  inelastic 
‘  scattering  Is  treated  as 
I  absorption. 

At  this  time  this  Indicator 
must  be  zero. 

Next  3  columns. 


General 

Description 


I  Format  Specific  Input  Quantity 

^  I sour! 

I  If  this  indicator  «  0,  a  table 
j  will  be  required  in  category  6 
input.  This  table  is  the 
I  spatial  source  at  the  first 
lethargy  level,  k-1. 

If  this  indicator  >  0  an  E. 
source  (as  described  in  the  text) 
is  computed  with  its  origin  at 
[  the  spatial  point  x  «  0.0. 

1  Next  3  columns. 


ID  10, 

If  this  indicator  ■  0,  there  is 
no  arbitrary  input  of  D^q  sources. 
If  this  indicator  0  an  Input 
table  FIST,  category  7,  will 
load  D,q  sources  directly  from 
input  for  each  lethargy  level 
from  KBGIHS  to  KIST. 

Next  3  columns. 


I  IROCO. 

I  If  this  indicator  ■  0,  the  compu- 
i  tation  of  the  flux  array  will 
f  proceed  by  row  iterations  only. 

{  If  IROGO  ■  0,  the  flux  computa- 
I  tion  will  proceed  alternately 
}  by  row  and  Column  iterations. 

1  Next  3  columns. 


ISYMM 

If  this  indicator  ■«  1  the  flux 
configuration  will  be  computed 
symmetric  about  the  left  edge 
of  the  array  (usually  x  ■  0). 

This  condition  of  symmetry  applies 
only  for  IROCO  «  0,  i.e. 
row  iteration  only  used  in  the 
computation  of  flux. 

Next  3  columns. 
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i 


Category  Mo. 


General 

Description 


Floating 

Point 

Input  Data 

Each  Single 
entry 
occupies  8 
columns. 

Each  input 
entry  of 
the  form 
.fff  X  10±®® 
is  expressed 
(.  ffffiH^e) 
in  the^S 
column  field. 

(A  certain 
econooy  is 
permitted  by 
Fortran  in 
preparing  E 
format 
Floating 
Point  input 
numbers. 

See  ITM 
Reference 
Manual,  709/ 
7090  FORTRAN 
PROGRAMMING 
SYSTEM  P.45, 
bottom. ) 


Format  Specific  Input  Quantity 


9E(8.2)  ALPLAY  (I),  I  -  1,  MIFLAG. 

This  table  defines  the  spatial 
and  lethargy  coordinates,  for 
the  finest  mesh,  of  the  flux 
array.  The  space  defining 
data  is  Ax^,  X2,  AX2 . 

*INONU*  ^INONU*  *IN0NU^^* 


"INONU*  ^INONU*  INONU^^* 

X  is  the  coordinate  of  the 
m 

lefthand  edge  of  the  mth  physi¬ 
cal  region;  Ax_  is  the  constant 

Ill 

interval  irlthin  that  region. 
Spatial  data  is  terminated  by  a 
constant  Interval  dlstance^O. 0. 
The  mesh  point  distances  give 
the  points  for  the  finest 
spatial  mesh:  hence  there  are 
J2  such  points.  The  table  con¬ 
tinues  with  the  beginning  entry 
of  the  lethargy  coordinates 
(usually  0.0).  Then  follows  the 
constant  Au,  the  maximum 
lethargy,  and  the  indicator 
value  0.0.  With  1  spatial 
region,  therefore,  there  would 
be  8  entries  or  64  columns  to 
this  table.  Where  INONU  no. 
of  regions,  use  next  16  INONU  + 
48  columns. 


criterion  for  flux  convergence. 
The  maximum  residual  obtained 
from  the  comparison  of  each  old 
flux  to  the  newly  computed  flux 
must  be  less  than  this  value  for 
convergence  to  occtir.  Next  8 
columns. 


is  a  very  small  positive 

nund^er.  The  number  of  flux 
values,  whose  absolute  value  in 
any  given  iteration  is  smaller 
than  is  cotmted  and  stored 

in  NCON  (location  53151g).  Next 
8  columns. _ 

ENZRI 

The  energy  at  which  lethargy 
li  «  0  for  the  particular  problem. 

Navt*  ortliifimo' 


51 


General 

Description 


Specific  Input  Quantity 


I  If  FLI  >  0,  all  values  of  the 

I  flux  taFle  are  initialized  to 
I  the  value  FLI. 

!  If  FLI  <  0,  the  flux  arrav^  \  f  \ 
!  k  is  initialized  to  e"l^k^ 

5  FLt  is  usually  set  equal  to 
‘  zero.  Next  8  columns. 


W  is  an  acceleration  factor. 
If  .-i  is  the  flux  computed 
on  tne'‘~(i+l)  St  Iteration, 
^i+1  replaced  by(W0j^_j_j^ 

+  (l-w)Cl.).  If  w«l,  this 
procedure  gives  nothing  new 
and  the  relationship  is 
skipped.  Next  8  columns. 


i  If  XO  >  0  the  left  edge 
5  boundary  of  the  flux  array  is 
I  initialized  to  the  value  of 


.! 


nuclide  2413 

identifica¬ 
tion  for 
region 


!  XO.  If  X  <  0  the  left,  or 
1  X  ,  boundary  fluxes  are  read 

i  in  by  input.  Next  8  columns. 

} 

j  XN. 

j  If  XN  >  0,  the  right  edge 
j  boundary  of  the  flux  array 
I  is  initialized  to  the  value 
i  of  Xjj.  If  Xjj  <  0  the  right 
;  hand  boundary  fluxes  are  read 
I  in  as  input.  Next  8  columns. 

f- - 

j  lELEW  (I)  I  -  1,24 
1  The  nuclide  identification 
numbers  for  all  nuclides 
within  the  first  spatial  region 
are  read  in  on  a  single  card, 
in  numerical  order.  These 
numbers  correspond  to  nuclide 
data  stored  on  tape  NUTAPE, 
the  nuclide  data  storage 
tape.  Each  entry  occupies  3 
columns.  This  table  takes  up 
n3  columns  of  one  card,  when 
n>number  of  nuclides  within* 
the  region  of  conqputatlon 


Category  No. 


General 

Description 


density  of 
nuclides 
for  a  given 
region 


format 


(E8.2) 


IDEBUG  is  a 
category  2 
indication 


Specific  Input  Quantity 


DENSIT(I), 

the  density  (gms/cc)  for  each 
iof  the  nuclides  above. 
iThis  pattern  of  a  single  card 
jof  lELEW  data  (4a  category 
I  card)  followed  by  the  corres- 
i ponding  density  card(s) 

I  (4b  category  card(s))ls 
!  maintained  until  all  INONU 
regions  of  the  problem  have 
their  composition  defined. 
Each  entry  has  8  columns . 

This  4b  table  occupies  a 
i  total  of  n8  columns,  where 
n^number  of  nuclides  within 
the  region  being  computed. 


optional  | 
input .  I 
IDEBUG-20  I 
indicates  I 
there  is  no  ! 
nuclide  data  | 
storage  tape  ! 
and  we  therefore 
read  category  f 
5  input .  I 


[  (E13.7, 

I  16,11) 


!  ENZR, 

I  the  energy  corresponding  to 
I  zero  lethargy  for  the  lethargy 
I  argument  table  of  the  cross- 
I  section  data.  This  must  be 
'  the  maximum  energy  for  the 
I  cross  sectiion  table.  First  13 
;  columns  floating  point  form. 
NOXSAG, 

I  the  number  of  entries  in  the 
lethargy  table  corresponding 
to  cross  section  input.  Next  6 
columns,  fixed  point. 

IBUGl, 

an  indicator.  Next  column  only. 
If  IBUGl-1  we  have  constant 
cross  sections.  Therefore 
we  read  aQ,aj^,  02,0^,73,7^  and 
73  tables  directly. 
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Category  No. 


General 

Description 


Format 


Specific  Input  Quantity 


Each  of  categories  5b,  5c>  and  5e 
through  51,  require  a  table  of 
numbers  to  be  prepared  for  input. 
Each  table  entry  is  in  floating 
point  form  where'  a  number  in  the 
form  +  .ffffff  X  10+®®  is  ex¬ 
pressed  +  .ffffffE+ee.  It  is 
placed  to  the  extreme  right  of 
the  13  column  field  reserved 
for  it.  One  IBM  input  card  can 
hold  up  to  4  entries.  (65  columns) 


15(13.7)  XSLEAG  (I)I-  1,  NOXSAG, 

the  table  of  lethargy  arguments 
:  for  all  cross  sections.  See 
note  following  5a^just  above. 


5c  i  Optional 
!  input;  if 
1  IBUGl-1, 

I  we  have 
I  constant 
I  cross  sec- 
I  tion  input, 
i  Therefore 
I  read  Input 
'  5  .  This 
I  option  is 
valid  only 
I  if  there 
{  is  1  region. 


i(El3.7)  !  The  following  will  form  tables 

;  with  constant  values 

I  XSMATS(l)  the  total  macroscopic 


scop; 

-13) 


cross  section  columns  (1 


ALFAO  (1)  value  of  a^^columns  (14-26) 
ALFA  1(1)  value  of  Oj^columns (27-39) 
ALFA  2(1)  value  of  a2Columns (40-52) 
ALFA  3(1)  value  of  columns (53-65) 
ALFA  4(1)  valu.  of 
GAM  3(1)  value  of  columns (14-26) 

GAM  4(1)  value  of  columns (27-39) 

GAM  5(1)  value  of  7^  columns (40-52) 

XSEE  91)  value  of  ^10^^^ 

columns (53-65) 


Category  No. 

General 

Description 

Format 

Specific  Input  Quantity 

So 

Optional 

Input.  If 
IBUGl  ^  1, 
we  have  non¬ 
constant 
cross  sec¬ 
tion  input 

given  in  5d- 
5e.  This 
option  per¬ 
mits  arbitrary; 
introduction 
of  SPM  I 

coefficients  ; 
without  refer- ; 
ence  to  NUTAPEj 
tape .  I 


5d 

i 

i 

\ - 1 

1  I7,3A6, 

2E15.6 

[  ■  1 
i 

i  . 

j 

f,  i 

l  ; 

1 

\ 

;  ■  ■  ■  ■  -  '  —  '  —  . . 

i  NUCLNO, 

'  the  identification  number  of 
the  nuclide  for  which  we  are 
reading  in  cross  section  data 
columns  (1-7) 

NAME,  a  18  character  title  for 
,  this  nuclide.  Columns  (8-25) 

AW,  the  atomic  weight  for 
this  nuclide,  columns  (26-40) 

1  USQ,  the  quantity  <p.^> 

:  columns  (41-65) 

5e 

j 

i 

1 

1  5E13.7 

I 

i  i 

{  i 

!  1 

•.  ! 

;  i 

:  \ 

1  XSMIT(I)  I»l,  NOXSAG 

1  o^,  or  microscopic  total  cross 
f  sections  in  barns 
;  See  note  following  category  5a 
i  description. 

5f  ; 

1 

1 

} 

/ 

1  •  i 

^  5E13.7  j 

i 

i 

! 

!  XSMIS  (I)  1=1,  NOXSAG 
j  a  ,  or  microscopic  scattering 

1  cross  section  in 

1  See  note  following  5a 

1  description. 

i 

5g  ; 

i 

\ 

\ 

1 

1 

(5E13.7)  1 

i  ; 

CIO (I)  I-l,  NOXSAG 

SPM  coefficient  defined  in  text. 
See  note  following  5a  descrip¬ 
tion. 

5h  ; 

1 

\ 

\ 

(5E13.7) 
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C20(I)  I-l,  NOXSAG 

Se®  note  following  5a  description 

Category  No. 


51 


5j 


5k 

51 


6 


7 


8 


9 


General 

Description 

Format 

Specific  Input  Quantity 

(5E13.7) 

C02(I)  I  -  1,  NOXSAG 

See  note  following  5a  description 

(5E13.7) 

Ci2(I)  1=1,  NOXSAG 

See  note  following  5a  description 

(5E13.7) 

WVO(I)  1=1,  NOXSAG 

See  note  following  5a  description 

' 

(5E13.7) 

WVl(I)  I  =  1,  NOXSAG 

See  note  following  5a  description 

Note  that  If  IDEBUG  ^  20,  all 
category  5  and  6  data  are  read 
from  the  nuclide  data  storage 
tape  NUTAPE. 

Optional 

Input . 

If  ISOUR  1 
=  0,  read 
the  lethargy 
=  0  source 

(9E8.2) 

FLIIX(I)  I  =  1,  J2 

Fluxes  at  lethargy  zero  level. 
Including  corner  points. 

J2  X  8  coliunns ,  with  at  most 

9  entries  on  a  card. 

Optional 

Input.  If 

DOlO  =  1, 
read  In 
tables  of 
arbitrary 
sources 

(9E8.2) 

FIST (I) 

Arbitrary  Dj^q  sources.  For 

each  lethargy  level  from  KBGIHS 
(category  2)  to  KIST  (category 

2),  Inclusive,  J2  (category  2) 
values  of  j  ^  will  be  read 

In,  from  x^  to  Xjj. 

(KIST-KBGIHS+1)  sets  of 
(J2  X  8  coltmms) 

Optional 

Input. 

If  XO  <  0 
read  In  X. 
side  ° 

boundary 
condition 

(9E8. 2) 

FLUX  (I) 

K2  entries  of  the  left  or  X 
boundary  values  of  flux, 
starting  from  lethargy  level  K=l. 
(K2  X  8)  columns 

Optional 

Input 

If  X„  <  0, 
read"ln 

XN  side 

boundary 

conditions 

(9E8.2) 

FLUX (I) 

K2  entries  of  the  rl^t  or  X^ 
boundary  values  of  flux,  ” 

starting  from  lethargy  level  K=l. 
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V*  SPM  Output 


The  SPM  output  of  results  Is  Initiated  by  a  printout  of 
a)  the  names.of  nuclides,  and  b)  their  corresponding  densities 
In  nuclei/cm-^  in  each  successive  region. 

By  a  proper  choice  of  the  Input  indicators,  IPCY  and  IMALP 
(category  2  Input) ,  the  SPM  program  may  be  set  up  to  print  out  con¬ 
vergence  data,  for  each  or  any  flux  Iteration  number.  In  the  following 
form; 

Line  1:  (flux  Iteration  number),  (number  of  points  tested  for 
convergence) ,  (average  residual) ,  (maximum  residual) , 

(j  position  of  max.  residual  (space)),  (k  position 
of  max.  residual  (lethargy)) 

Line  2:  (Bl,  Minimum  Value  of  Flux),  (B^  Maximum  Value  of  Flux), 
(absolute  value  of  maximum  of  (B2,-Bl)) 

At  iteration  IMALP  and  at  completion  of  the  problem  by  either  con¬ 
vergence  or  divergence,  SPM  prints  out  problem  identification 
comments.  SPM  briefly  prints  out  why  the  problem  ended.  It  then 
prints  the  flux  array  requested  by  input  indicators  MKO  and  MJO 
(category  2  input) .  The  energy  level  and  space  position  identifying 
each  element  of  the  flux  array  is  printed  on  each  page. 


VI.  GENUC  -  Cross  Section  Processing  Program 

A.  Purpose 

The  purpose  of  the  GENUC  program  Is  to  produce  NUTAPE 
from  microscopic  cross  section  data.  It  produces  also  an  edited 
version  of  NUTAPE  for  printing.  A  printed  record  Is  necessary 
for  the  operator  to  have  In  order  to  Identify  the  various 
nuclides.  The  following  data  are  given  on  NUTAPE  for  each 
nuclide.  A,  Og,  C^q,  C2Q,  0^2*  Cj^2*  '^o*  '^l* 

has  an  option  whereby  Cj^2  ''^1  equal  to  zero.  This 

Is  equivalent  to  dropping  the  third-order  terms.  This  option  Is 
denoted  In  the  following  by  7  =  0. 

All  data  thus  collected  for  a  given  nuclide,  along  with 
title  Information,  are  written  on  tape  as  a  single  record.  In  the 
exact  order  SPM  requires. 

GENUC  provides  for  the  sequential  storing  of  many 
nucllde-data  records  on  the  same  tape.  A  necessary  restriction 
of  storage  on  one  tape  is  that  each  nuclide  must  use  the  same 
table  of  lethargy  arguments.  The  lethargy  table  and  title 
Information  to  Identify  this  particular  tape  are  found  as  the 
first  record. 

B.  Tape  Unit  Designations 

GENUC  consists  of  a  main  program  and  one  subroutine 
CHANTP.  Because  of  the  organization  of  the  subroutine  CHANTP, 
the  GENUC  code,  like  the  SPM  program,  may  be  run  at  any  7090 
Installation  regardless  of  the  standard  Input  and  output  tape 
designations  of  that  Installation. 
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The  GENUC  code  deck  Issued  with  this  report  uses  the 
logical  tape  units  listed  In  the  table  below. 


Unit 


Use 


2  Input 

3  Output,  Hollerith 

5  Output  of  nucllde-data,  In  binary  form, 

for  use  with  the  SPM  program 

The  procedvire  for  changing  the  logical  tape  designations  that 
GENUC  uses  Is  exactly  the  same  as  for  the  SPM  code. 


C.  Input  Preparation  and  Restrictions 

An  Input  deck  for  the  preparation  of  a  new  data¬ 
storage  tape  (l.e.  when  the  Input  quantity  N00NTP  >=  0)  must 
carry  both  a  lethargy  division  table  and  Its  associated  base 
energy  for  u  >=  0,  both  of  which  will  be  maintained  as  a  standard 
for  all  successive  lethargy-dependent  variables  on  this  tape. 

When  N00NTP  >  0,  the  lethargy  table  will  be  assumed  present 
within  the  first  tape  record  and  the  Input  deck  need  only  con¬ 
tain  information  directly  related  to  the  given  nuclide.  Any 
number,  N0NUC,  of  different  nuclides  may  be  added  to  the  tape, 
with  two  restrictions;  1)  GENUC  provides  no  interpolation 
facilities  for  non-standard  lethargy  arguments;  and  2)  the  use 
of  the  7“  option  (defining  7  to  be  zero  or  non- zero)  remains  the 
same  for  all  N0NUC  nuclides.  When  the  GENUC  run  Is  complete  for 
these  N0NUC  records,  the  GENUC  program  starts  anew.  At  this  time 
for  the  new  run,  the  7  option  of  course  may  be  changed. 
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GENUC  DEFINITIONS 


N0NUC 

N00NTP 

IBUG 

XSLEAG(I) 

N0XSAG 

NIDTAP 

ENZR 

NUCLN0 

NAME 

AW 

USQ 

FL(KK,L) 

XSMIT(L) 

XSMIS(L) 


Number  (Index  maximum)  of  nuclides  to  be  added  onto 
data  storage  tape 

The  actual  number  of  nuclide  records  on  the  nticlide 
data  tape.  If  0,  then  we  prepare  a  new  nuclide  tape. 

Indicator  for  7  option 

When  IBUG  =1,  the  tables  C12(L)  and  WVl(L)  are  set 
equal  to  zero  for  all  lethargies  L. 

Table  of  lethargy  arguments  for  cross-section  tables 
with  a  N0XSAG  number  of  entries. 

The  number  of  entries  to  above  table,  <  200. 

A  label,  or  title,  to  be  printed  out  for  identification 
of  a  particular  nuclide  data  tape 

The  base  energy  for  which  lethargy  is  defined  as  zero 
Number- label  of  a  nuclide 

18  Alpha-numeric  character  (3  cell)  title  for  a  nuclide 
Atomic  weight  of  the  nuclide 


Legendre  coefficients  f.  for  the  nuclide.  Here 
KK  =  i  +  1  ^ 

There  are  7  separate  such  tables  for  KK  =  2  to  KK  =  8 
L  is  the  lethargy  index 

Each  such  table  has  N0XSAG  entries  in  it,  corresponding 
to  the  XSLEAG  table  of  arguments 

Microscopic  cross-section  in  barns  corresponding  to 
XSLEAG  lethargy  table,  N0XSAG  entries 

Microscopic  cross-section  in  barns  corresponding  to 
XSLEAG  lethargy  argument  table,  N0XSAG  entries 
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OUTPUT  OF  GENUC  PROGRAM 


The  First  Record  has 
NIDTAP 
N0XSA6 

XSLEAG(I),  I  •=  1,  N0XSAG 

Each  successive  record  describes  a  distinct  nuclide  and  contains: 
NUCLN0 

MAME  18  characters 
AH  atomic  weight 
XSMIT(L)  o^ 

XSM1S(L)  a  scattering 

C10(L)  (Note  that  defined  In  SPM  text,  Is 

C20(L)  never  distinct  from 

C02(L) 

*  C12(L) 

VO(L) 

*  VI  (L) 
end  of  record 


*  We  note  that  where  IBUG  -\(G1  input),  the  C12(L)  and  also  the 
VI (L)  tables  will  be  set  equal  to  zero. 
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